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Abstract 

We present the first unquenched, continuum limit, lattice QCD results for the matrix el- 
ements of the operators describing neutral kaon oscillations in extensions of the Standard 
Model. Owing to the accuracy of our calculation on AS = 2 weak Hamiltonian matrix 
elements, we are able to provide a refined Unitarity Triangle analysis improving the bounds 
coming from model independent constraints on New Physics. In our non-perturbative com- 
putation we use a combination of Nf = 2 maximally twisted sea quarks and Osterwaldcr- 
Seiler valence quarks in order to achieve both 0(a)-improvement and continuum- like renor- 
malization properties for the relevant four-fermion operators. The calculation of the renor- 
malization constants has been performed non-perturbatively in the RI-MOM scheme. Based 
on simulations at three values of the lattice spacing and a number of quark masses we have ex- 
trapolated/interpolated our results to the continuum limit and physical light/strange quark 
masses. 
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1 Introduction 



The fundamental target of present-day research activity in Particle Physics is the search for 
New Physics (NP) effects beyond the Standard Model (SM) predictions. Two main routes 
are followed, one based on the study of processes in which the direct production of NP 
particles at high energy colliders like the LHC takes place, and a second one based on the 
indirect investigation of NP effects coming from the exchange of virtual NP particles. 

In the so-called indirect approach a crucial role is played by Flavor Physics processes 
that are sensitive to NP through loop effects. These processes vanish at tree level in the SM, 
and some of them are theoretically very clean, despite the fact that they are loop mediated, 
and in some cases also CKM or helicity suppressed. Among them, AF = 2 transitions have 
always provided some of the most stringent constraints on NP. For instance, the constraints 
from K° — K° oscillations are particularly stringent for NP models that generate transitions 
between quarks of different chiralities [1]-[4 |. Therefore, an accurate determination of the 
AS = 2 bag parameters (S-parameters) is crucial to the improvement of NP constraints. 

In the present work, we provide the first accurate lattice determination of the AS = 2 
i?-parameters relevant for physics beyond the SM, calculated in the continuum limit and 
using data from unquenched, Nf = 2, dynamical quark simulation^. Our results represent 
a significant improvement with respect to the (quenched) input values of Refs. (6| and Q 
used so far in phenomenological analyses. The result, obtained within the same lattice 
setup, for the Bk parameter, relevant for K° — K° mixing in the SM, has been presented 
in®. 

The outline of the paper is as follows. Section [2] contains a brief description of the 
AS = 2 matrix elements of the effective weak Hamiltonian describing the most general 
pattern of K° — K° oscillations. In Section [3j based on the results of this work for the 
AS = 2 B-parameters, we discuss the implications for NP of our updated Unitarity Trian- 
gle (UT) analysis Q. In Section |4] we illustrate the main theoretical features of the lattice 
setup employed in our simulations (twisted mass lattice QCD [10], iTTTlD and we describe 
the strategy for obtaining accurate numerical estimates of the ^-parameters as well as ra- 
tios of kaon four-fermion matrix elements. In Section [5] we collect our numerical results. In 
Section|6]we give our estimates of the various ^-parameters and matrix elements ratios and 
compare the present results with the previous determinations existing in the literature. In 
five Appendices we discuss a number of technicalities: i) the renormalization properties of 
the four-fermion operators in our "mixed action" setup lfl2~l : ii-iii) the RI-MOM computa- 
tion of renormalization constants (RCs) and corresponding results respectively; iv) tables of 
lattice data on pseudoscalar meson masses, decay constants and bare four-fermion matrix 

2 While finalizing our paper we became aware of the work of Ref. where B-parameters and ratios of 
four-fermion matrix elements have been computed at one lattice spacing with Nf — 2+1 dynamical quarks. 
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elements; v) complete results on renormalized four-fermion matrix elements obtained by 
using various formulae for the chiral extrapolation and two alternative procedures for the 
RI-MOM determination of RCs. 

2 AS = 2 effective weak Hamiltonian 

The general form of the AS = 2 effective weak Hamiltonian is 

5 3 

^e A / =2 = E Ci<Di + C % 6 t , (2. 1) 

i=l i=l 

where in the so-called SUSY basis (12, |[T3l ) the four-fermion operators Oi and Oi have 
the form 

O x = [^ 7M (l- 75 )d Q! ][^(l-75)^] 

2 = [ s -(l- 75 )^][^(l- 75 )^] 

3 = [i«(l- 75 )^][^(l- 75 )d a ] (2.2) 

04 = [s a (l- 7 5)^][^(l+ 75 )^] 

5 = [^(1-75)^(1+75)^] 

01 = [s a 7M (l + 75 )<i Q ][s /3 7, t (l + 75)^] 

02 = [s a (l + 75)^^(1 + 75)^] (2-3) 

03 = [^(1 + 75)^] [^(1 + 75)^] 

with a and /3 denoting color indices. Spin indices are implicitly contracted within square 
brackets. The Wilson coefficients Cj and Ci have an implicit renormalization scale depen- 
dence which is compensated by the scale dependence of the renormalization constants of 
the corresponding operators. 

Notice that the parity-even parts of the operators Oi and Oi are equal. From now on and 
for notational simplicity we will denote by Oi (i = 1, . . . , 5) the parity-even components 
of the operators (12.21) . Due to parity conservation in strong interactions, in the study of 
K° — K° oscillations it is then sufficient to consider only the matrix elements (K°\Oi\K °). 
We recall that in the SM only the kaon matrix element of the operator 0\ comes into play. 

The bag parameters, Bi (i = 1, . . . , 5), provide the value of four-fermion matrix ele- 
ments in units of the magnitude of their vacuum saturation approximation. More explicitly 
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they are defined by the equations lfl4l 



(K^O^K ) = ^Bxi^mlfl (2.4) 
(K^O^K ) = ZiBiQjt) \ m *l K - J 2 fari = 2 5, (2.5) 

with£j = (8/3, —5/3, 1/3, 2, 2/3). In the above relations one recognizes B\ as the famil- 
iar Bk- We recall that, as suggested by the parametrization adopted in the r.h.s. of Eqs. (12.41 ) 
and (12.5I ). the (K \Oi (/x) | K° ) matrix element is expected to vanish in the chiral limit unlike 
the other four ones. 

An alternative way which has the merit of allowing a more accurate evaluation of 
matrix elements, is to consider the matrix elements ratios Ri = (K°\Oi\K ) / [K Q \0\ \K°), 

1 = 2, . . . , 5, as first proposed in Ref. For details on our lattice implementation see 
Section @] and in particular Eq. ( 14.19I ). 

For the reader's convenience we here anticipate our final continuum results for Bi and 
Ri in the MS scheme of Buras et al, defined in Ref. iTTSl . and the RI-MOM scheme at 

2 GeV, see Tables Q] and ^respectively. Details on the calculation and uncertainty estimates 
are given in Section [5] In Appendix |E] we also give the final continuum results for Bi and 
Ri in the MS and the RI-MOM scheme at 3 GeV. 



MS (2 GeV) 



Bi 


B 2 


B 3 


Ba 


B 5 


0.52(02) 


0.54(03) 


0.94(08) 


0.82(05) 


0.63(07) 




Ri 


R2 


Ri 




Rs 


l 


-14.7(06) 


5.2(04) 


25.7(11) 


6.8(05) 



Table 1: Continuum limit results for Bi and Ri, renormalized in the MS scheme of Ref. lfl5l 
at 2 GeV. 



3 Model-independent constraints on AS = 2 operators and New Physics scale 
from the Unitarity Triangle analysis 

AF = 2 processes provide some of the most stringent constraints on NP generalizations of 
the SM. Several phenomenological analyses of AF = 2 processes have been performed in 
the last years, both for specific NP models and in model-independent frameworks. A gener- 
alization of the UT analysis, which allows for NP effects by including the most significant 
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RI-MOM (2 GeV) 



Bi 


B 2 


B 3 


Ba 


B 5 


0.51(02) 


0.73(04) 


1.29(11) 


1.04(07) 


0.76(09) 




Ri 


R2 


Ri 




Rs 


l 


-13.6(06) 


4.9(03) 


22.5(10) 


5.7(04) 



Table 2: Continuum limit results for Bi and Ri, renormalized in the RI-MOM scheme at 
2 GeV. 



flavour constraints on NP available at the time was performed in Ref. [9]. The result was 
a simultaneous determination of the CKM parameters and the size of NP contributions to 
AF = 2 processes in the neutral kaon and B^ s meson sectors. 

The NP generalization of the UT analysis consists in including in the theoretical parametriza- 
tion of the various observables the matrix elements of operators which, though absent in the 
SM, may appear in some of its extensions. The analysis shows that the constraints coming 
from K° — K° matrix elements are the most stringent ones, in particular for models that 
generate transitions between quarks of different chiralities (see Refs. [fl]] - (4l). Thus an 
accurate determination of the AS = 2 S-paiameters is crucial to the improvement of the 
NP constraints. 

The results for the AS = 2 £>-parameters obtained in the present work come from 
unquenched Nf = 2 lattice QCD data carefully extrapolated to the continuum limit. They 
hence represent a significant progress with respect to the input values used in the UT analy- 
sis performed in Ref. [0, where quenched lattice numbers (without a systematic continuum 
limit extrapolation analysis) computed more than five years ago in Refs. [ 6 ] and were 
employed. For this reason, we present here an update of the analysis of Ref. [9] based on 
our new values of the AS = 2 i?-parameters. The new ingredients entering the analysis 
are collected in Tables Q] and |2l For all the other input data we use the numbers quoted in 

Ref. m. 

In the present NP-oriented analysis, the relations among experimental observables and 
the CKM matrix elements are extended by taking into consideration the most general form 
of the AS = 2 effective weak Hamilonian (see Eq. (12.11) ). The effective weak Hamilonian 
is parameterized by Wilson coefficients of the form 

Q(A) = ^, i = 2,...,5, (3.1) 

where F{ is the (generally complex) relevant NP flavor coupling, Lj is a (loop) factor which 
depends on the interactions that generate Cj(A), and A is the scale of NP, i.e. the typi- 
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cal mass of new particles mediating AS = 2 transitions. For a generic strongly interact- 
ing theory with an unconstrained flavor structure, one expects Fj ~ Lj ~ 1, so that the 
phenomenologically allowed range for each of the Wilson coefficients can be immediately 
translated into a lower bound on A. Specific assumptions on the flavor structure of NP 
correspond to special choices of the Fi functions. For example Minimal Flavor Violation 
(MFV) models El-Ell correspond to F 1 = F SM and F^ = 0. 

Following Ref. [9], in deriving the lower bounds on the NP scale A, we assume Li = 1, 
that corresponds to strongly-interacting and/or tree-level coupled NP. Two other interesting 
possibilities are given by loop-mediated NP contributions proportional to either a 2 s or a w . 
The first case corresponds for example to gluino exchange in the minimal supersymmetric 
SM. The second case applies to all models with SM-like loop-mediated weak interactions. 
To obtain the lower bound on A entailed by loop-mediated contributions, one simply has to 
multiply the bounds we quote in the following by a s (A) ~ 0.1 or aw ~ 0.03. 

In agreement with Ref. [9], we find that in the K° sector, due to the non- vanishing 
chiral limit (chiral enhancement) of their matrix elements, all bounds coming from the con- 
tributions of non-standard operators (i.e. from the operators 0{ with i ^ 1) are more than 
one order of magnitude stronger than the bound from the SM 0\ operator. 





95% allowed range 


Lower limit on A 




(GeV- 






(TeV) 


ImCf 


[-2.8,2.6] 


• io- 


15 


1.9 


• 10 4 


ImCf 


[-1.6,1.8] 


• 10" 


17 


24 


10 4 


Im Cf 


[-6.7,5.9] 


• io- 


17 


12 


10 4 


ImCf 


[-4.1,3.6] 


• io- 


18 


49 


10 4 


Im Cf 


[-1.2,1.1] 


• io- 


17 


29 


10 4 


ImCf 


[-4.4,2.8] 


• io- 


-15 


1.5 


• 10 4 


ImCf 


[-5.1,9.3] 


• io- 


17 


10 


10 4 


Im Cf 


[-3.1,1.7] 


• io- 


16 


5.7 


• 10 4 


ImCf 


[-1.8,0.9] 


• io- 


17 


24 


10 4 


Im Cf 


[-5.2,2.8] 


• 10- 


17 


14 


10 4 



Table 3: 95% probability range for the lmC l K coefficients and the corresponding lower 
bounds on the NP scale, A, for a generic strongly interacting NP with generic flavor structure 
(Li = Fi = 1). In the lower panel the results of (9l are displayed for comparison. 

The results for the upper bounds on the ImCf coefficients and the corresponding 
lower bounds on the NP scale A are collected in Table [3] where they are compared to the 
previous results of Ref. [9]. The superscript K is to recall that we are reporting the bounds 
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50 




Ci C 2 C 3 C 4 C 5 



Figure 1 : The lower bounds on the NP scale, provided by the constraints on Im Cf (i = 
1 , . . . , 5) for generic NP flavor structure, are shown as brown bars. For comparison, we plot 
the bounds of Ref. [9] as yellow bars. 

coming from the kaon sector we are here analyzing. Although several input parameters 
have been updated with respect to Ref. @ (see Ref. [16]), the more stringent constraints 
on the Wilson coefficients of the non-standard operators and, consequently, on the NP 
scale, mainly come from the improved accuracy achieved in the values of the AS = 2 
^-parameters obtained in the present work. This can be realized by comparing the small 
improvement of the bound coming from Im Cf , obtained using a value of the B\ -parameter 
very close to the one taken in Ref. (9J, with those coming from the other coefficients using 
the new i?-parameters. We observe that the analysis is performed (as in (9]) by switching on 
one coefficient at the time in each sector, thus excluding the possibility of having accidental 
cancellations among the contributions of different operators. Therefore, the reader should 
keep in mind that the bounds may be weakened if, instead, some accidental cancellation 
occurs. 

In Fig. [T] we show the comparison between the lower bounds on the NP scale obtained 
for the case of a generic strongly interacting NP with generic flavor structure by the con- 
straints on the Im Cf coefficients coming from the present generalized UT analysis, and 
the previous results of Ref. (9) . 

As a specific example of NP models we consider the warped five-dimensional exten- 
sions of the SM discussed in Ref. ll23l . where the origin of hierarchies in quark masses 
and mixings is explained via the localization properties of quark wave functions in the fifth 
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dimension. In particular, in the Randall-Sundrum (RS) scenario one has 

U = (g* s f, F 4 = ^f^, A = Mq , (3.2) 

where Mg and g* ~ 6 are the mass and coupling of Kaluza- Klein excitations of the gluon, 
Y* ~ 3 is the five-dimensional Yukawa coupling (whose flavour structure is assumed to be 
anarchic), ~ 3 MeV and m s ~ 50 MeV are MS quark masses at the high scale and 
v = 246 GeV is the Higgs vev. Running from a reference scale of 5 TeV, we obtain at 95% 
probability \mCf G [—6.5, 8.2] • 10~ 18 , from which we get 

M G >49TeV. (3.3) 

Considering instead gauge-Higgs unification (GHU) models, from Ref. [23] we have 

L, = (g* s f, F^^f^, A = Mq , (3.4) 

where in this case g* ~ 4 is the five-dimensional gauge coupling in units of the radius of 
the compact dimension. We obtain the bound 

M G >74TeV. (3.5) 



4 Non-perturbative lattice computation of the AS = 2 matrix elements 

Lattice QCD provides an ideal first principle framework in which non-perturbative compu- 
tation of hadronic matrix elements can be performed with controlled systematic uncertain- 
ties. In the last few years a number of lattice determinations of Bk of increasing precision 
have appeared in the literature based on a variety of lattice regularizations with Nf = 2 or 
N f = 3 dynamical fermions OD, EU-HH. For recent reviews see Refs. l29l-|[3T1l. 

Very little has been done in the literature concerning the calculation of the physical 
matrix elements of the full operator basis of Eq. (12.2b - Calculations of the whole set of 
AS = 2 renormalized operators have been carried out, in the quenched approximation, 
using improved Wilson fermions (see Refs. [6] and [ 14]) or the chirality conserving overlap 
and domain-wall fermions (see Refs. (7J and 11331 . respectively). Very recently a N* = 2 + 1 
dynamical quark calculation appeared Q that uses domain-wall fermions at one value of 
the lattice spacing^]. 

The first calculation of the kaon matrix elements of the whole operator basis was per- 
formed employing Clover improved Wilson fermions. Since the Clover term coefficient was 

3 A preliminary computation of the bare matrix elements using unquenched Nf — 2 + 1 domain-wall 
dynamical fermions was presented in Ref. 1321 . 
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set to its tree-level value, matrix elements were affected by 0(<?Qa) discretization errors. 
Simulations were carried out at two values of the gauge coupling corresponding to lattice 
spacings ~ 0.07 and ~ 0.09 fm (0, lfl4l0 . The major source of systematic errors was, how- 
ever, the uncertainty related to the construction of the multiplicatively renormalizable lattice 
operators Oj. In fact, owing to the breaking of the chiral symmetry intrinsic in the Wilson 
fermion action, the bare counterparts of each of them mix with all the other operators of 
equal dimension including those with "wrong" chiral transformation properties Il34l . All 
the mixing coefficients and the overall RC were computed in the non-perturbative RI-MOM 
scheme ll35l . 

Similar quenched computations were carried out using overlap and domain-wall fermions 
(see Refs. [7] and [33], respectively). Though performed at pretty coarse lattice spacings 
(namely a ~ 0.09 and ~ 0.13 fm in the case of overlap fermions and a ~ 0.1 fm in the 
case of domain-wall fermions), these simulations have the advantage that the renormaliza- 
tion properties of the operators entering the four-fermion basis are as in the continuum, and 
lattice artifacts are 0(a 2 ). Also in this case the non-perturbative RI-MOM scheme was used 
in the computation of the various RCs. 

In the following sections we present a new unquenched computation of the kaon ma- 
trix elements of the full AS = 2 four-fermion operator basis employing lattice data from 
simulations with Nf = 2 dynamical fermions, performed at three rather fine values of the 
lattice spacing in the interval [0.07, 0. 1] fm. We are thus able to safely extrapolate the lattice 
estimators of all the relevant matrix elements to the continuum limit (CL). 

We use a mixed fermion action setup where we adopt different regularizations for sea 
and valence quarks. In particular we introduce maximally twisted (Mtm) sea quarks iTTTTl 
that we take in combination with Osterwalder-Seiler (OS) [36] valence quarks. This strat- 
egy has been suggested in Ref. |[T2l as a way of setting up a computational framework 
allowing for a calculation of AS = 2 four-fermion matrix elements that is both automati- 
cally O(a) improved and free of wrong chirality mixing effects. A proof of the latter point 
is given in Appendix lAl and its validity is numerically verified in Appendix [Bj while O(a) 
improvement of physical quantities is a genuine property of the setup of Ref. lfl2l . As a 
consequence unitarity violations due to different sea and valence quark regularization yield 
only 0(a 2 ) artifacts, provided renormalized sea and valence quark masses are matched. 
In our case, the matching of the renormalized quark masses is obtained by simply taking 
identical values for the corresponding sea and valence bare mass parameters. 

The interesting lattice setup briefly described above has already been successfully 
tested in Bk computations both in the quenched approximation [37 ] and on ensembles 
with Nf = 2 [8] and Nf = 2 + 1 + 1 dynamical quarks II38T1 . as well as in unquenched 
(Nf = 2 + 1 + 1) studies of meson masses and decay constants [39] and nucleon sigma 
terms |30). 
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4.1 Sea and valence quark regularization 



The Mtm-LQCD action of the light quark flavor doublet can be written in the so-called 
"physical basis" in the form ifTTI 

ST" = « 4 £ \ £ 7m(V m + V* ) -i 75 r 3 r sea [M CI - 1 £ V m V m] . 

(4.1) 

The subscript sea is to remind us that this action will be used to generate unquenched 
gauge configurations. The field tp describes a mass degenerate up and down doublet with 
bare (twisted) mass fi sea - The parameter M cr is the critical mass that one has to fix non- 
perturbatively at its optimal value (as proposed in Refs. HTTl - 11431 and implemented in 
Refs. [44] and [45 ]) to guarantee the 0(a)-improvement of physical observables and get 
rid of all the unwanted leading chirally enhanced cutoff effects. In the gauge sector the 
tree-level improved action proposed in Ref . [46 ] has been used. 

For valence quarks we use the OS regularization [36]. The full valence action is given 
by the sum of the contributions of each individual valence flavour qj and reads Ifl2l 

S? a S i = « 4 E ^ X ){\ E >(Vm+V* )-vysr f [M CI ~ £ V^V M ] +^}q f {x) , (4.2) 

where the index / labels the valence flavors and M cr is the same critical mass parameter 
which appears in Eq. (14- 1 b - We denote by tf and fj,f the values of the Wilson parameter and 
the twisted quark mass of each valence flavor. 



4.2 Lattice operators and correlation functions 

In the strategy proposed in Ref. lfT2l . which we follow here, four species of OS valence 
quark flavors (qt, f = 1, . . . , 4) are introduced, two of which (qi and q%) will represent the 
valence strange quark with masses fii = fi^ = /U« s », while the other two (q2 and q^) will be 
identified with the light up/down quarks having masses fi2 = = fii- The corresponding 
r f Wilson parameters must obey the relation 

r\ = r 2 = r 3 = -r 4 . (4.3) 

In the numerical computations reported in the present work we have averaged over the two 
cases r\ = ±1, holding r 2 , and related to r\ as in Eq. (14.3b - 

As in the case of the computation of -Ba-(= B\), in the calculation of the bag parame- 
ters Bi (i = 2, . . . , 5), we need to consider the axial currents 

A f = <7i7m75<72 Af = (?37^75<?4- (4.4) 
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and the pseudoscalar quark densities 



P LZ = <?i75<72 P U = <7375<?4 

In addition, we need to consider the following set of four-fermion operators 

Ofi = 2{([ff7M9S1[^7Mflf] + [$1^MW<&\) ± (2 o 4)} 



(4.5) 



4 f4 = 2{([g?g 2 Q ][g^] " [#769?] 



^75<fi)±(2o4)} 



MA 

5[±] 



2{([3f«J][8ftf] " [^76^][^76fl?]) ±(2^4)} 



M)±(2^4)} 
-f 75 gf])±(2o4)} 

3' 



(4.6) 



where square parentheses denote spin invariants and a and /3 are color indices. 

In the mixed action (MA) approach defined above the following properties can be 
proved (see Appendix A and Ref. lfl2l ). 

(i) The operators 0¥^ defined in Eq. (14.61 ) enjoy continuum-like renormalization prop- 
erties, 



I u i[+] \ 

U 2[+] 
U 3[+] 

n MA 

\ U 5[+] I 

where the matrix is defined in Eqs. (IA.7l ). (IA.9I ) and (IA. 10b . 



/ Z n 














\ 


( U A+] 


\ 





Z22 


•^23 










nMA 

U 2[+] 







Z32 


•^33 










nMA 

U 3[+] 













Z44 


Z45 




nMA 

U 4[+] 




V o 








Z54 


Z55 


) 


V U 5[+] 


I 



(4.7) 



(ii) The axial currents and pseudoscalar quark densities, defined in Eqs. (14.41 ) and (14- 5b . 
are renormalized according to the formulae ( fl2l . fiTTl ) 



{ pl2 



121 

fi Jren 



Z A A 
Z S P 



12 

12 



■"■y. Jren 
.34] 



[P 



ZvAf. 
Z P P U 



(4.8) 
(4.9) 



(iii) The matrix elements (P 43 |0^f |P 12 ) reil , obtained from correlation functions of the 
renormalized operators in Eqs. (I4.7I )- (|4.9I ), with the identification \i\ = n$ = /j, s (= 
bare strange quark mass) and [12 = jt/4 = (= bare up-down quark mass), will tend 
in the limit a — > to the continuum matrix elements (K°\Oi\K°) of the (parity-even 
parts of the) operators of Eq. (12.21 ) with mere 0(a 2 ) discretization errors. 
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Ref. GO 



In the construction of correlation functions we follow the general procedure outlined in 
(H. We use periodic boundary conditions in every direction for all fields, except for the 
quark fields on which we impose anti-periodic boundary conditions in the time direction. At 
time slices yo an d yo+T/2 "wall" operators with K°-meson quantum numbers are inserted. 
The first operator is constructed in terms of q% and q\ quark fields and the second in terms 
of ^4 and qs quark fields. Explicitly they have the expressions 

y 

Vy* + T = (|) 3 ^g4(y,yo + r/2)7 5(? 3(y,yo + T/2) (4.10) 

y 

The correlators we then need to compute are0 

Ci(x ) = (|) 3 E^'o-ff °f4^' x °) V vl) > < = 1, ■ ■ ■ , 5 , (4.1 1) 

Cpp(x ) = (|) 3 £(P 12 (x,*o)0> (4.12) 

X 

C' PP (x ) = (f) 3 E^ 4 o+f pU & x °))- ( 4 - 13 ) 

X 

To improve the signal-to-noise ratio a sum has been performed over the spatial position of 
the four-fermion operator, and for each gauge configuration the time slice yo is randomly 
chosen. An important contribution to the reduction of statistical fluctuations comes also 
from summing over the spatial position of both kaon interpolating fields in Eq. (14.111) . After 
summing over the the spatial position of the four-fermion operator, which is known from 
experience to be crucial for the signal, to project on zero 3-momentum states, just one 
further spatial sum is needed. The second spatial sum, which we are able to do, gives a 
further signal improvement. These spatial sums were implemented and carried out at a 
reasonably low computational price by means of the stochastic technique discussed in sect. 
2.2 of Ref. QQ. 

For large time separation yo xq <C yo + T/2 the (plateau of the) following ratio 
estimate 

E ^ 1w - c^M^r i = 2 -" 5 (414) 



For the special case of i = 1, the evaluation of Bi = Bk (see Ref. 1 8 1) requires the use of the two-point cor- 
relation functions that involve the axial current and have the form Ca p (xo) = (a /Lj X^(^o 2 (x, x o) "Pyl ) 

mdC' PAo (x ) = (a/L) 3 E 3 <^ 4 3 + x A?(3,xo)). 
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provides an estimate of the (i = 2, . . . , 5) bag parameter since 

y «x «y +T/2 (K°\0^\K°) {b) 



(K°\P 12 \0) (0|P 34 |K°) 



i t Bf\ i = 2 ) ...,5, 

(4.15) 



with &, U, ft) = (-5/3, 1/3, 2, 2/3). 
We also compute the ratios 



E[R? ] \(xo 



Cj(x ) 
Ci(* )' 



2,. ..,5. 



(4.16) 



which for large time separations yield the ratios of the corresponding kaon four-fermion 
matrix elements, according to the formula 



(K°\0^\K°) 



(6) 



R 



(b) 
i ' 



2,. ..,5. (4.17) 



4.3 Estimates of renormalized quantities 

Recalling Eqs. (14.71) and ( 14.91 ), the renormalized values of the bag parameters will be given 
by the formula 

where a sum over the repeated index j is understood. As for the renormalized expression of 
the four-fermion operator ratios of Eq. (14. 17b . we choose to evaluate the rescaled quantity 

- _ ,f K y r M 12 M 34 z tJ (6)] 

where M /fc and -F ifc are the mass and decay constant of the pseudoscalar meson made of 
the (qiqk) quark pair. In order to compensate the chiral vanishing of the {K°\Oi\K°) 
matrix element we have multiplied the ratios R.f^ by the factor M 12 M 34 /F 12 F 34 ; its form 
has been chosen in such a way to partially compensate the lattice artifacts affecting the 
different lattice discretizations of kaon mesons (resulting from different choices of the OS 
r^-parameters) we use. Furthermore, we have remultiplied the quantity in the square bracket 
by the ratio of the experimental values of the kaon decay constant (/|? p = 156.1 MeV) and 
its mass (m^ p = 494.4 MeV). The definition we get in this way is in line with the one 
proposed in Ref. [7 ]. Based on the discussion of Section 4.2 (in particular item iii), we note 
that in the continuum limit and at the physical values of the u/d and s quark masses the 



5 In the following the superscript (b) denotes bare quantities. 
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quantity R{ of Eq. (14.191 ) provides the right estimate for the ratio of the renormalized matrix 
elements of interest, i.e. 

Rl ~ (K°\ Ol \K0) (4 - 20) 

We end this section by recalling that in the twisted mass mixed action setup of Ref . lTT2l 
we are using, lattice estimators of physical quantities are only affected by 0(a 2 ) lattice 
artifacts. This is true in particular for the kaon mass and decay constant as well as for the 
kaon four-fermion matrix elements. 



5 Simulations, data analysis and results 

The ETM Collaboration has generated Nj = 2 configuration ensembles at four values of 
the inverse bare gauge coupling, j3, three of which will be used in the present work, and 
at a number of light quark masses, /x sea . The values of the simulated lattice spacings lie 
in the interval [0.07, 0.1] fm. Bare quark mass parameters are chosen so as to have light 
pseudoscalar mesons ("pions") in the range 280 < mpg < 500 MeV and heavy-light 
pseudoscalar mesons ("kaons") in the range 450 < mps < 650 MeV. Simulation details 
are given in Table [4] 

The value of the light u/d quark mass parameter, afj,g, is common to sea and valence 
quarks, while the heavier quark (the would-be strange quark that we denote by "s", see 
Table |U) is quenched. As discussed in Section 15.21 we will get to the physical kaon mass 
by suitably interpolating (extrapolating) data in fj,« s » (m) to the "physical" value /j, s (/J, u /d), 
while simultaneously taking the continuum limit. The "physical" values fi u u and [a s of 
the quark masses are known and can be found in Ref. 11481 . The quark bilinear RCs, Zp 
and Zg, have been computed in the non-perturbative RI-MOM scheme in Ref. H71 . A RC 
computation for the full basis of the four-fermion operators using RI-MOM techniques is 
presented in Appendix IB1 In Appendix ICl we collect the values of the four-fermion RCs that 
are used in this work, as well as Zp and Z$- 



5.1 Extracting bare estimates from lattice data 

Bare results for the ratio of the four-fermion matrix elements (c.f. Eq. (14.171 )) with i = 
2, . . . , 5 at the three (3 values and combinations of quark masses are listed in Tables [141 [T31 
and [16] of Appendix iDl In Tables [T7l[T8l and [191 of the same Appendix we collect the results 
for the bare quantities & Bf (i = 2, . . . , 5) at each value of f3. 

For illustration in Fig. [2] we display some examples of the B f plateau quality at f3 = 
3.8, 3.9, 4.05 for the smallest value of the light quark mass (panels (a), (d) and (c)), and 
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/3 = 3.80, a~0.10fm 


a^i£ — ap, sea 


a" 4 (L 3 x T) 


a/i« s >' 


Ngtat 


0.0080 


24 3 x 48 


0.0165, 0.0200, 0.0250 


170 


0.0110 


it 




180 


P = 3.90, a ~ 0.09 fm 


0.0040 


24 3 x 48 


0.0150, 0.0220, 0.0270 


400 


0.0064 


it 




200 


0.0085 


a 




200 


0.0100 


a 




160 


0.0030 


32 3 x 64 




300 


0.0040 


66 




160 


p = 4.05, a ~ 0.07 fm 


0.0030 


32 3 x 64 


0.0120, 0.0150, 0.0180 


190 


0.0060 


it 


a 


150 


0.0080 


it 




220 



Table 4: Details of simulation runs at P = 3.80, 3.90 and 4.05. 



at P = 3.90 for the next-to-smallest quark mass (panel (b)). Vertical dotted lines indicate 
the plateau region where the K°- and the T^-state dominate the three-point correlators. 
Similar examples of the plateau quality for the case of the four-fermion operator ratios are 
illustrated in Fig. [3] Both Figs.[2]and[3]display very good signals. 

5.2 Computation at the physical point 

Extracting physical quantities from lattice data requires performing extrapolations and/or 
interpolations of renormalized lattice estimators to the physical point (continuum limit and 
"physical" value of quark masses). 

5.2.1 RCs computation and combined continuum- chiral extrapolation 

We have computed the full matrix of the four-fermion operator RCs in a mass independent 
scheme. We carry out the non-perturbative calculation adopting the RI-MOM approach. 
The implementation of the RI-MOM setup has been presented in Refs. [8] and (47). We 
should mention that in our RC estimators cutoff effects, though parametrically of 0(a 2 ), are 
numerically reduced owing to the subtraction of perturbatively evaluated 0(a 2 g 2 ) contribu- 
tions. After that, two different, but by now standard fl47l . procedures are employed to deal 
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Figure 2: Data and time-plateaux for E[Bf >}(t) (i = 2,..., 5) (see Eq. KIM ) plot- 
ted vs. 2t/T = 2(xo — yo)/T. In panel (a) we show data for (3 = 3.80 and 
(ane,afi V ) = (0.0080,0.0165) on a 24 3 x 48 lattice; in panel (b) for (3 = 3.90, 
(ane,afi V ) = (0.0040,0.0150) on a 24 3 x 48 lattice; in panel (c) for (3 = 4.05 and 
{ani,ati« a ») = (0.0030,0.0120) on a 32 3 x 64 lattice; in panel (d) for j3 = 3.90 and 
(afj,£,afj," s ") = (0.0030,0.0150) on a 32 3 x 64 lattice. Vertical dotted lines delimit the 
plateau region. For clarity data for E[B^](t) have been slightly shifted. 



with 0(a 2 p 2 ) discretization effects. The first, called Ml, consists in linearly extrapolating 
to zero the residual (after the perturbative subtraction) 0(a 2 p 2 ) terms. The second one (so- 
called p 2 -window method, or M2 for short) leads to RC estimates obtained by averaging 
data over a fixed (in physical units) and very narrow momentum interval. 

We carry out continuum and chiral extrapolations in a combined way. For all bag 
parameters, Bi, and ratios, R{ (see Eq. ( 14.191) ), we have tried out a fit ansatz of the following 
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Figure 3: Data and time-plateaux for E[R\ (see Eq. S4.16\) ) plotted vs. 2t/T = 
2(xq — yo)/T; panels (a), . . ., (d) correspond to cases with i = 2, . . . , 5 respectively. In 
each panel we show data for /3 = 3.80 and (a/^, a^< V ' ) = (0.0080, 0.0165) on a 24 3 x 48 
lattice; (3 = 3.90, (a^,ayu« s ") = (0.0040,0.0150) on a 24 3 x 48 lattice; ft = 4.05 and 
(a/^,a/i« s ») = (0.0030,0.0120) on a 32 2 x 64 lattice. Vertical dotted lines delimit the 
plateau region. 



general form 

2 

Y = ^ 4 n W s ) [r„ hT + D Y {r Q jx s ) [-f , (5.1) 

n=0 r ° 

where we have made explicit the dependence of the fit parameters Ay ' and Dy on the 
renormalized strange quark mas^f] in units of vq (ro/x s ). We studied separately the cases 
of linear and polynomial ansatz. We have also considered NLO ChPT fit functions for Bi 

6 We use the symbol (~) to denote renormalized quark masses in the MS scheme at 2 GeV. 
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based on the formulae given in [49] in the case of SU(3). Those formulae transformed to 
NLO SU(2) ChPT read: 

+ D Bi (r ^)\-} 2 (5.2) 

with B = 2.84(11) GeV (renormalized in MS at 2 GeV) and /„ = 121.0(1) MeV, as we 
used in [8]. The sign before the logarithmic term is minus (-) for i = 1, 2, 3 and plus (+) for 
i = 4,5. As for Ri and i = 2, 3 the ChPT fit formula at NLO coincides with the linear fit 
ansatz, while for the cases i = 4, 5 we use 

+ D Ri {r^ s )[yJ (5.3) 

The (ro / a) values at each (3 are 

(Wa)/3=3.80 = 4.54(07), (r /a)p =3 . 90 = 5.35(04), (r /a)^ =4 . 5 = 6.71(04) (5.4) 

The u/d and s quark masses have been computed in Ref. H§1 . Their values in the MS 
scheme at 2 GeV are 

/uJ|(2GeV) = 3.6(2) MeV, ^(2 GeV) = 95(6) MeV (5.5) 

In the four panels of Fig. [4] we show the combined chiral and continuum fit (see 
Eq. (15.11 ) and Eq. (15.31 )) for the ratios Ri against the renormalized light quark mass for 
i = 2, . . . , 5, respectively. The RCs used in these plots are the ones computed with the Mi- 
method and are expressed in the MS scheme of Ref. |[T5l at 2 GeV. Lattice data correspond 
to points taken at the pair of quark masses (rofte, rop, s ). 

In panels (a) and (b) (corresponding to cases with i = 2, 3 respectively) we display 
the curves that correspond to the polynomial fit function (15.11 ) at the three j3 values we are 
considering in this paper. The black solid line represents the continuum limit curve. The 
dashed black line represents the continuum limit curve that is obtained if a linear fit ansatz in 
fi£ is used. Black open circles and triangles stand for the results at the physical quark mass 
point from the polynomial and the linear fit ansatz, respectively. Recall that for Ri with 
i = 2, 3 ChPT fit formula at NLO coincides with a linear fit ansatz. In panels (c) and (d) 
(corresponding to cases with i = 4, 5 respectively) we display the curves corresponding to 
the ChPT fit formula and the linear fit function. In this case black open circles and triangles 
stand for the results at the physical quark mass point from the ChPT fit and the linear fit 
ansatz, respectively. 

Similarly, in Fig. [5] we present the combined chiral and continuum fit for the Bi- 
parameters, again renormalized in the MS scheme of Ref. IIT51 at 2 GeV. In all four panels 



Bi = B*(roA«) 



1 + bi(r fi s ) =F 



^log 



2(4vr/ ) 



(4tt/ ) 



Ri = Rf(r fi s 



1 + Cj(roAs) + 



^log 



(4tt/, 



oj 



(4tt/ q ) 5 
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Figure 4: Solid lines in panels (a) and (b) show the behaviour vs. the renormalized light 
quark mass of the combined chiral and continuum fits (according to the polynomial for- 
mula H5. iD with n = 2) of the Ri ratios, with i = 2 and i = 3 respectively, renormalized 
in the MS scheme of Ref. £/J/ at 2 GeV with the Ml -type RCs. In panels (c) and (d), 
solid lines, instead, show the combined chiral and continuum fit described by NLO-ChPT, 
Eq. (T5.3D for i = 4 and i = 5, respectively. The full black line is the continuum limit 
curve. The dashed black line represents the continuum limit curve in the case of the linear 
fit ansatz. Black open circles and triangles stand for the results at the physical point corre- 
sponding to the polynomial (panels (a) and (b)) and ChPT fit (panels (c) and (d)), and linear 
fit ansatz, respectively. 



we display the curves corresponding to the ChPT fit formula at the three f3 values. The 
black solid line represents the continuum limit curve. The dashed black line represents the 
continuum limit curve that is obtained if a linear fit ansatz in p,? is used. 
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Figure 5: Solid lines in panels (a) to (d) show the behaviour vs. the renormalized light 
quark mass of the combined chiral and continuum fits (according to the ChPT fit formula 
(5.2\) ) for the Bi parameters with i = 2, . . . , 5 respectively, renormalized in the MS scheme 
of Ref. 051? at 2 GeV with the Ml-type RCs. The full black line is the continuum limit 
curve (f5. JD . The dashed black line represents the continuum limit curve in the case of the 
linear fit ansatz. Black open circles and triangles stand for the results at the physical point 
corresponding to the ChPT fit and linear fit ansatz, respectively 



Note the nice agreement (within one standard deviation) of the two fit ansatze for both 
the Ri ratios and the Bi parameters. 

In Tables [5] and [6] we gather our final continuum resulta^ for Ri and Bi in the MS of 



7 They are the same results presented also, for reader's convenience, in Tables [JJ and [2] in the introductory 
section[2] 
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Ref. HH and RI-MOM scheme at 2 GeV. The final value of B t for i = 2, . . . , 5 has been 
computed by averaging the estimates obtained from the three kinds of fit ansatz discussed 
above, and using bootstrap error analysis. The half difference between the two more distant 
results has been taken as an estimate of the systematic error associated to the extrapolation 
procedure. The total uncertainty is obtained by adding in quadrature the statistical and the 
systematic error. For i = 1 we present the result for Bk published in Ref. JSJ. 

In Appendix |E] we provide more detailed results obtained from the various fitting pro- 
cedure we have investigated. We also show that the continuum extrapolated quantities that 
are eventually obtained by employing Ml-type and M2-type RCs turn out to be perfectly 
consistent between each other within statistical errors, except for one case where cut-off 
effects on M2-RCs are too big for a significant comparison. Also in Appendix |EJ see Ta- 
bles [22] and [23} we quote our continuum results for Bi and R, L in the MS and the RI-MOM 
scheme respectively at 3 GeV. 

As already stated, an alternative (indirect) way to compute the ratio of the kaon matrix 
elements of the renormalized operators Oj, i = 2, . . . , 5 to that of 0\ is based on the formula 
(see eqs. (I2~4b and (|23T » 

(K°\Oi<ji)\K°) eiUi(/i) (Ml*) + /^)) 2 ' 

This of course requires knowledge of the Bi parameters and the renormalized quark masses. 

We find that the two evaluations (indirect and direct, based on Eq. (15.61 ) and Eq. (14. 19b . 
respectively) lead to compatible results within errors. However, the indirect estimates suffer 
from much larger final uncertainties. This is due to several reasons. One is related to the 
quadratic dependence on the quark mass, which makes the relative error on the mass to give 
a significant contribution to the final error. Furthermore in the indirect method one has to 
consider extra uncertainties due to the error of the bilinear operators' RCs that are used to 
compute the £?-parameters. A comparison of the direct and indirect results obtained for the 
ratios i?j is provided in Appendix |E] 
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MS (2 GeV) 



Bi 


B 2 


B 3 


Ba 


B 5 


0.52(02) 


0.54(03) 


0.94(08) 


0.82(05) 


0.63(07) 




Ri 


R2 


Ri 




Rs 


l 


-14.7(06) 


5.2(04) 


25.7(11) 


6.8(05) 



Table 5: Continuum limit results for Bi and Ri, renormalized in the MS scheme of Ref. Ifl5l 
at 2 GeV. 



RI-MOM (2 GeV) 



Bi 


B 2 


B 3 


B A 


B 5 


0.51(02) 


0.73(04) 


1.29(11) 


1.04(07) 


0.76(09) 




Ri 


R2 


Rs 




R 5 


1 


-13.6(06) 


4.9(03) 


22.5(10) 


5.7(04) 



Table 6: Continuum limit results for Bi and Ri, renormalized in the RI-MOM scheme at 
2 GeV. 

6 Conclusions 

Accurate measurements of the K° — K° mixing amplitudes can yield useful hints on New 
Physics if theory can provide comparatively accurate calculations of quantities parametriz- 
ing beyond the SM effects. This requires a precise, first principle, evaluation of the kaon 
matrix elements of the full basis of four-fermion operators entering the most general effec- 
tive AS = 2 weak Hamiltonian. 

In this paper we have presented the first unquenched lattice QCD determination in the 
continuum limit of the matrix elements of the full AS = 2 four-fermion operator basis. 
We have used Nj = 2 unquenched tm-LQCD gauge configurations produced by the ETM 
Collaboration in combination with maximally twisted valence quarks of the OS type. 

The mixed action setup proposed in Ref. lfT2l offers the possibility of obtaining au- 
tomatically O(a) improved results and an operator renormalization pattern identical to that 
of a chirally invariant regularization at the rather cheap price of mere 0(a 2 ) unitarity vi- 
olations. Using data at three lattice spacings (with a in the interval [0.07, 0.1] fm) and a 
number of pseudoscalar masses ("pions") in the range [280, 500] MeV, we are able to safely 
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carry out the continuum and the light quark mass limit of the observables of interest. All 
results are non-perturbatively renormalized in the RI/MOM scheme. 

We get in this way the most accurate estimates to date of AS = 2 effective weak 
Hamiltonian matrix elements. The total error on the Ri ratios is between 4% and 7% and 
on the bag parameters, B{, between 4% and 12%. 

Tables |7] and [8] show a comparison between our results for Ri and B{ (in RI/MOM at 
2 GeV) and the data at fixed lattice spacings coming from the two old quenched calculations 
of Refs. El and OlB 

For the ^-parameters, one finds large differences between the central values of our re- 
sults and those of Refs. [6 ] and Q, which vary between 10% and 30% (though the errors are 
typically comparably large). With respect to Ref. [6], the differences are even larger when 
the results are compared in terms of the ratios Ri, presumably due to a combined effect, 
in this case, of having overestimated the values for both B\ and the strange quark mass in 
the computation of [6]. We emphasize that, with respect to the old quenched calculations, 
having performed in the present study simulations at three values of the lattice spacing and 
quite smaller values of the pion masses provides us with a much better control over the main 
sources of systematic uncertainties, besides the quenched approximation. Current experi- 
ence suggests that the possible systematic errors related to the quenching of the strange and 
charm quarks, which still affect our calculation, are negligible within the present uncertain- 
ties. Forthcoming results from simulations in the continuum limit with Nf = 2 + 1 and 
Nf = 2 + 1 + 1 dynamical flavours will provide a check of this expectation. We should add 
that our continuum limit results for Rj and Bi (i = 2, . . . , 5) are in the same ballpark with 
the numbers given at one lattice spacing in Ref. (H where Nf = 2 + 1 dynamical quarks 
are employed. 

As an interesting phenomenological application of the results obtained in this paper 
we have carried out a new UT analysis along the lines of the work of Ref. Thanks 
to the improved accuracy of the present determination of the AS" = 2 5-parameters, we 
could substantially strengthen the existing upper bounds on the Wilson coefficients of the 
operators of the non-standard sector of the effective weak Hamiltonian, and consequently 
increase the lower bound on the New Physics scale. 



8 In this comparison we do not include the (preliminary) quenched results at one value of the lattice spacing 
given in Ref. 1 3 3 1. 
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Bi (RI-MOM at 2 GeV) 


This work 


Ref. Q 


Ref. El 


CL 


a = 0.09 fm 


a = 0.13 fm 


a = 0.07 fm 


a = 0.09 fm 




1 


0.51(02) 


0.56(05) 


0.53(04) 


0.68(21) 


0.70(15) 


2 


0.73(04) 


0.87(07) 


0.90(10) 


0.67(07) 


0.72(09) 


3 


1.29(11) 


1.41(12) 


1.53(40) 


0.95(15) 


1.21(10) 


4 


1.04(07) 


0.94(05) 


0.90(13) 


1.00(09) 


1.15(05) 


5 


0.76(09) 


0.62(05) 


0.56(14) 


0.66(11) 


0.88(06) 



Table 7: Comparison between the unquenched results for Bi obtained in the present work 
and the quenched values of Refs. Q and |6|. Numbers are for renormalized quantities in 
the RI-MOM scheme at 2 GeV. 



Ri (RI-MOM at 2 GeV) 


This work 


Ref. Q 


Ref. H 


CL 


a = 0.09 fm 


a = 0.13 fm 


a = 0.07 fm 


a = 0.09 fm 




1 


1 


1 


1 


1 


1 


2 


-13.6(06) 


-16.1(30) 


-15.8(29) 


-6.7(18) 


-6.6(11) 


3 


4.9(03) 


5.2(09) 


5.4(08) 


1.9(05) 


2.3(04) 


4 


22.5(10) 


20.7(30) 


18.8(28) 


12.1(33) 


12.6(21) 


5 


5.7(04) 


4.6(06) 


3.9(13) 


2.6(07) 


3.3(05) 



Table 8: Same as in Table [7] for the Ri ratios. 
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A Renormalization properties of AS* = 2 four-fermion operators 



In this appendix we want to spell out the renormalization properties of the four-fermion 
operators of interest for the description of K° — K° oscillations in the mixed action (MA) 
lattice setup of Section|4] We will do this by exploiting the results of Ref. [50]. In particular, 
we show that the operators in Eq. (14.6b . that we report here for the reader convenience, 



exhibit the same renormalization pattern as the corresponding continuum operators. We 
recall that the Wilson r-parameters of valence quarks in Eq. (IA. lb are taken as specified in 
Eq. (14.31) . The normalization we have chosen in the definitions (1A.1I) is such that, when the 
operators are taken between the pseudoscalar operators P 12 = (7175(72 and P 43 = 

<?475<?3> o ne gets the same Wick contraction multiplicities one would obtain in QCD upon 
evaluating the kaon matrix elements of the operators (12.2b - Naturally, apart from the issue 
of renormalization, the physical matrix elements will be obtained (in the continuum limit) 
by finally setting in our MA setup [i\ = /J3 = fj, s and [12 = /M = A% with fi s (resp. \X() 
corresponding to the bare strange (resp. degenerate up-down) quark mass. 

A key result of Ref. ll50l (stated there in the quark basis that is most natural for un- 
twisted Wilson fermions) is that the combinations that enjoy simple renormalization prop- 
erties are not the operators listed in Eqs. (IA. lb . but those that are obtained after performing 
a Fierz transformation on the operators 0^_^ and O^fy . This transformation has the effect 
of rewriting O^A and 0^_^ in terms of operators where both spin and color indices are 
contracted within the same pair of quarks. With the definitions 




Hiltfl^M"/^} + [tfWsqZMWscfi]) ± (2 o 4)} 



2{([M[«] + [$1^2Ml<\) ± (2 o 4)} 

2{([5?rf][gftf] + \%h*<g\[4'i*tf\) ± (2 o 4)} 

" Wls<M][4lwi\) ± (2 O 4)} 

i{{w4M<&] - w?is4Mis<&\) ± (2 o 4)} , 



(A.l) 




2{([<?i 7^2] [937^94] + [qi 7^ 7592] [937^75^4]) ± (2 o 4)} 

2{ ([<Zi7m<?2] [937^94] - [^17^75^2] [937^75^4]) ± (2 o 4)} 



mivm^qsqi] ~ [917592] fe7594]) ± (2 4)} 

^{{[qmWqm] + [qnm]%imi]) ±(2^4)} 

2{([9i<V92][930/ii/g4]) ± (2 O 4)} (for/i > u), 



(A.2) 
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where = [7^,7^/2, one gets 



n MA - \i ± 1n MA 



i[±] 



A [±] 



/ 1 














1 



Tl/2 



1 













±1/2 





(A.3) 



V Tl/2 

In Eq. (IA.2I ) we have omitted color indices as they are always contracted within each square 
parenthesis. 

In order to make direct contact with the formulae of Ref. [50] we must pass from the 
gj-basis, in which the valence quark action (14.21 ) was written and where the Wilson term is 
(maximally) twisted, to the x/-basis, where the Wilson term takes its standard form. This 
is achieved by the chiral transformation 



If 



-iirr /75/4 



9/, Qf 



Xf = 9/e 



-i7rr/75/4 



under which 

Q OS — V ci* os 

'-'val ' ^val 



(A.4) 



(A.5) 



x < - 



2 ^ 



7m(V m + V*) + [M cr - I V m V m] + *76»-/A*/}x/( 



and assuming r4 = ±1 (as used in the present work), 

Qn±] = 2 H(Rl7MX2][X37A t 75X4] + [Xl7^75X2][X37MX4]) ±(2h4)} 
= 2i {([WX2][X37M75X4] - [Xi7/,75X2][X37mX4]) t(2^4)} 

= 2i{([xi75X2]foX4] - [X1X2][X375X4]) T (2 o 4)} 
Qf[i] = 2i{([xi75X2]foX4] + [X1X2][X375X4]) ± (2 O 4)} 
2i{([xi^X2][X3^75X4]) ± (2 o 4)} . 



r.O MA 
r 4Vl[±] 

r 4V2[±] 
r,<7 MA 

r,<7 MA 

r,(7 MA 



f\MA 
V5[±] 



(A.6) 



According to Ref. 11501 . for the renormalized operators Qf /A one g ets 



I ^i[±] 

nMA 
^2[±] 

nM4 

^3[±] 

^4[±] 
nMA 



> ren 


/ Z11 














\ 


[±] 




\ 







^22 


^23 












r,MA 

( *2[±] 









^32 


-^33 












/)MA 
^3[±] 















Z44 


-^45 






nMA 

^4[±] 




/ 


V 








-2<54 


-£55 


) 




nMA 
\ ^5[±] 


/ 



(b) 



(A.7) 



Since at fi f = the fermion action is indistinguishable from a standard massless Wil- 
son fermion action, in any mass independent renormalization scheme, the operators Q 



MA 

i 
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in the l.h.s. of Eqs. (1A.6I) enjoy the same renormalization properties of the corresponding 

operators Qf IA into which they are trasformed under (1A.4I) . The Qf IA operators have, up 

to lattice artefacts, the same RCs of the corresponding operators in the standard Wilson's 

formulation of lattice QCD. This result could also have been proved using the somewhat 

more elaborated approach of Ref. Ifl2l . 

From Eqs. (IA.6MA.7I) and recalling Eq. (IA.3I) . one finally arrives for the operators 
i 

[+] 



Off,i of interest to us in this paper at the renormalization formulae 



U i[+] 



z = a [+] z q (aI 



Z ij U j[+] 



+h-i 



(6) 



(A.8) 
(A.9) 



Zr 



\ 




















^22 


-Z H 

^23 











^"32 


^33 



















7 {+] 

^"45 











crl+] 
^"54 


^"55 



/ 



(A. 10) 



At this point the matrix elements (P 43 |0^| rcn |P 12 ) | rcn , built using Eqs. dA.lt . dA.8t 
and (IA. 10b and evaluated in our MA setup with /ii = /X3 = (J, s (strange quark mass) and 
^2 = M4 = W (up-down quark mass), tend in the limit a — > to the matrix elements 
(K°\Oi\K°) of the operators dZ2t in QCD with mere 0(a 2 ) discretization errors ifTHl . 
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B RI/MOM computation of renormalization constants of four-fermion oper- 
ators 



In order to convert our lattice results for the bag parameters Bi to their physical continuum 
counterparts, in the same renormalization scheme and at the same scale as the correspond- 
ing perturbative Wilson coefficients used in the phenomenological analysis, we need the 
renormalization constants (RCs) of the operators 0¥A , i = 1, 2, . . . , 5, see Eq. (14.6b . or 
equivalently Eq. (IA.1I ). As discussed in Section 4, in our mixed action (MA) setup for lat- 
tice correlation functions these operators represent the analogs of the parity-even parts of 
the AS = 2 four-fermion operators (12.21 ) that are relevant in the formal continuum theory. 

In this Appendix, we give details on the non-perturbative computation of the RCs 
performed using the well established, and by now standard, RI-MOM [35] scheme 0. 

As explained in detail in Ref. ll50l . instead of using the operator basis 0¥_$ of Eq. (IA.1I ). 
it is more convenient to employ the Fierz transformed operators Q¥_^, i = 1, 2, . . . , 5, de- 
fined in Eq. (IA.2b . To lighten our notation, in the following we will drop the superscript and 
the sign subscript, denote these operators simply by Qi and assemble them in the array Q. 
The generic renormalization pattern of the bare operators is of the form 

Q rcn = Z [I + A] QW (B.l) 

where the scale-dependent renormalization matrix Z is block-diagonal, with a continuum- 
like block structure (the same as for e.g. the matrix in Eq. (1A.7I )). while A is a sparse 
off-diagonal and scale-independent matrix of the form 






A12 


A13 


A14 


A15 


A 2 i 








A24 


A25 


A 3 i 








A34 


A35 


A 4 i 


A42 


A43 








A51 


A52 


A53 









However, as shown in Appendix A, using the MA lattice setup of Section 4, the "wrong chi- 
rality mixing" terms Ajj, are reduced to mere 0(a 2 ) effects, the renormalization matrix Z 
coincides with the matrix Zq of Eq. (1A.10I ) and we recover a continuum-like renormaliza- 
tion pattern. This is a very important advantage of our approach, which we will implement 
in practice using the following strategy: compute the quark propagators in the gy-basis 
(also called physical basis of tmQCD at maximal twist, in which the critical Wilson term 
is twisted, see Eq. (14.21 )). impose RI-MOM renormalization conditions on the operators Qi 
and extract the renormalization matrix (Z) and, for check purposes, the mixing matrix (A). 

'Actually we use the RI'-MOM scheme. The prime on RI is to remind about the specific definition we 
adopted for the quark field RC, Z q = £1, where £1 is the quark propagator form factor defined in Eq. (3.5) of 
Ref. (47). 
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B.l Procedure for extracting the RCs 

To determine the matrices Z and A in Eq. (IB. lb we proceed as follows. We start by com- 
puting the lattice quark propagator 

S qf (p) = a A Yl e ~ tpX <?/(*) ?/(0)> (B.3) 

X 

and the four-point Green functions with an insertion of the operator Qi, namely 

Gi(p, p, p, pTJ^s = (B.4) 

a 16 £ e - lp(xi - X2+X3 - Xi) { [fi(xi)] . fe(x2)] 6 Q . (Q) [ft(xj|)] c [ft(a . 4)] {) . 

El ,112,2:3 ,^4 

The lower( upper) case Greek (Latin) symbols denote uncontracted spin (color) indices. The 
corresponding amputated Green functions are given by 

Ai(j>,P,P,p)£$ s = (B.5) 

For the sake of clarity, we will use matrix notation, denoting the matrices by boldface sym- 
bols and omitting color and spin indices. The amputated Green functions will be collected 
in the 1x5 row vector 

A(p) = (Ai, A 2 , A 3 , A 4 , A 5 ) (p,p,p,p) . (B.6) 

Setting 

k(ap,afj) = Z- 2 (ap)A(ap) [I + A T ] Z(o/i) T , (B.7) 

the renormalization matrix Z(a / u) is determined by solving the renormalization condi- 
tions 1501 . namely 

PA(p) | p2 = M 2= I. (B.8) 

In Eq. (IB.8I ). Z q is the quark field RC and A is, as we said before, the mixing matrix. We 
have also introduced the 5x1 column vector of spin projectors (see Eq. (37) of Ref. iTSUt 
for the explicit form of these projectors) 

P T = (Pi, P2, P3, P4, P5) (B-9) 

which act on the amputed Green functions by (i,j = 1 • • • , 5), 

PiAj = Tr PiAj(ap) 
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where the trace is taken over spin and colour, and obey the orthogonality relations 

TtPiAf (ap) = <%, (B.10) 

with the tree level amputated Green function of the operator Qj. It is convenient to 
express A in terms of a "dynamics" matrix D, defined by 

A(ap) = A (0) (ap) D(ap) . (B.ll) 

This matrix equation can be solved for D using the spin projectors P, getting 

D(p) = PA(p). (B.12) 

Combining Eqs. (IB.7I) . (IB.8I) and (IB.12I) . we see that, once the dynamics matrix is known, 
we can determine both the renormalization and the mixing matrices from the relation 

Z q 2 D [I + A T ] Z T = I — > Z [I + A] = Zg (D T ) _1 . (B.13) 

This matrix equation can be solved for Z and A by exploiting the block diagonal structure 
of the Z matrix. In fact, it is easy to see that the three diagonal blocks of the renormalization 
matrix (i, j = 1; i , j = 2, 3 and i, j = 4, 5) are given by 

Zij = Z 2 q ( D T )~ 1 = 1) (», j = 2, 3) (i, j = 4, 5) (B.14) 

whereas the mixing coefficients are easily obtained from the equations, 

i = 2,-.- ,5 (B.15) 

i = 2 j = 1,4,5 1 
i = 4 i = 1,2,3 J ' 

We can now summarize our procedure to determine the renormalization matrix of the 
parity-even part of the four-fermion operators of the SUSY basis of Eq. (12.21) . 

Step 1 The Green functions (IB .41 ) and (IB.5I) . and from them the dynamics matrix D, are 
evaluated in the Landau gauge for a sequence of sea, [i sea , and valence, \i va i, quark 
mass values at each of the three lattice spacings we consider here. The bare param- 
eters and the statistics of this computation are detailed in Table 2 of Ref. H71l . One 
can also find there (see Eqs. (3.6) and (3.7)) the set of discrete lattice momenta, p u 
(Pi,2,3 = (27r/L) 711,2,3, P4 = (27r/T) (714 + 1/2)), that we include in the present 



Zu A 



11 



Zi 



Hi+1 



A; 



TT 1 



A 



i+lj 



(D T ) 
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calculation. To minimize the contributions of Lorentz non-invariant discretization 
artifacts, we take into consideration only momenta satisfying the constraint 

X>p < 0.28 y52f v J, ap v = sm(ap„). (B.16) 

In the following, we shall often use the short-hand p 2 = pt- 

Step 2 For each f3 and each choice of the scale p 2 , the renormalization relation (IB. 141) is 
enforced at all values of \i sea and \i va \ given in Table 2 of Ref. 821 • By doing so, we 
obtain at nonzero quark masses the estimators Z^j 1 (p 2 ; a 2 p 2 ; fi va i; ^sea), which are 
then extrapolated to [i va i = (see Section IBT21 and fi sea = (see Section IB3T ). 

Step 3 Improved estimates of Z^ 1 ' (p 2 ; a 2 p 2 ;0; 0), are obtained by subtracting the pertur- 
batively leading cutoff effects (see Section uT4l . 

Step 4 Using the NLO continuum QCD evolution of the renormalization matrix Z calcu- 
lated in Refs. |[5Tl[T5l . the first argument of Zf^' (p 2 ;a 2 p 2 ; 0; 0) is brought to a refer- 
ence scale /ig. In this step, we assume that the scales p 2 and /ig are large enough to 
make NLO perturbation theory accurate. This is the same level of accuracy achieved 
in the determination of the Wilson coefficients. 

Step 5 The residual a 2 p 2 dependence in Zf^' (/Uq; a 2 p 2 ; 0; 0) is attributed to lattice artifacts, 
which we treat according to either the Ml or M2 methods, introduced in Ref. [47 ] (see 
Section IB. 61 ). 

Step 6 In order to reduce the statistical error, the lattice RC estimators are averaged over 
two equivalent patterns of Wilson parameters (n,r2, r^, r^), namely (1,1,1,-1) and 
(—1,-1,-1,1), as well as over different lattice momenta corresponding to the same 
p 2 . We have checked that performing these averages before or after taking the chiral 
limit leads to consistent results. 



B.2 Valence chiral limit 

In view of the relation (IB . 14b and since the extraction of Z g (see Ref. B71 ) poses no partic- 
ular problems, our discussion will be mainly focused here on the quark mass dependence of 
the dynamics matrix D. At fixed values of /3, a 2 p 2 and a^i sea , we fit the dynamics matrix 
elements Dij to the ansatz 

) = A(p 2 ;(i sea ) + B(p 2 ;fj, sea ) (i va i + C(p 2 ;n sea )/fj, va i. (B.17) 
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Here we have introduced a term with a pole in /j, va i ~ rnps to co P e w * tn tne ex P ect ed 
Goldstone boson (GB) pole contribution to the elements of the D matrix. The existence 
of such a GB-pole term can be understood as follows. At asymptotically large p 2 , non- 
perturbative effects giving contributions potentially divergent in the chiral limit to the Green 



functions (IB.4I) do vanish and the latter turn out to be polynomial in the quark mass parame- 
ters 051 . At finite values of p 2 , however, the contributions to (the spectral decomposition of) 
these Green functions from one-GB intermediate state with momentum q and mass mps, 
give rise to terms proportional to (q 2 + m'p s )^ 1 and suppressed by some power of 1/p 2 . 
If several one-GB intermediate states contribute to the spectral representation of the Green 
functions (IB .41 ) several terms, each behaving as {q 2 + m 2 PS )~ l and suppressed by some 
power of 1/p 2 , will show up. These results follow straightforwardly from the "polology' 
study of the Green functions (see e.g. the discussion in the book fl52"Tl ) or from the well 
known Lehmann-Symanzik-Zimmermann (LSZ) reduction formalism. Now, since in the 
Green functions (IB.4I ) the four-fermion operator is inserted at zero four-momentum transfer 
(q = 0), one expects, from the time orderings where two quark fields can create from the 
vacuum a pseudoscalar (i.e. GB) one-particle state, a contribution proportional to l/m 2 ps , 
suppressed by some power of 1/p 2 . Similarly, from those time orderings where two quark 
fields create a GB-state and two further quark fields destroy another GB-state, contributions 
do arise that behave as {l/m 2 PS ) 2 and are twice more strongly suppressed at large p 2 . 

In conclusion, by exploiting (along the lines of Appendix A of Ref. ll35l ) the large- 
p 2 behaviour of the matrix element [^1 (0\qf(p)qf/(—p)\pf'f) and taking also into consid- 
eration the four factors of S g (p)^ 1 that stem from the relations (IB.5I ). one finds that the 
dynamics matrix Dij contains GB-pole contributions of the following kinds: 

D Singlcl2 ^ P t S-^ p )S- 2 \-p)(0\q 2qi \P 12 } \ (P 12 \QmMPW-p)) 

p [M ) 

jingle 34 ^ Pl { qi (p)q 2 (-p)\Q ^P^) 1 (P^ \q m I05" 1 (p)S- 1 (~p) 

[M ) p 

n Doublcl234 pn-lyn-1/ n \/C\\nn I p!2\ (P 12 \Qj(fy\P U ) 

x (P 34 |g 4 g3|0)5- I (p)5- 1 (-P) • (B.18) 

We recall that the kinematics of our Green functions corresponds to an exceptional 
momentum configuration where p\ = p% = —p 2 = —p^ = p, and thus q = 0. As in 
the chiral limit S q ' 1 (p) ~ 7| t p^, the result (IB. 18b implies that single and double GB-pole 
terms are suppressed by 1/p 2 and l/(p 2 ) 2 factors, respectively. A second observation is 



10 Here q'f (p) (qf (p)) denotes the four-dimensional Fourier transform of the quark field q'f (x) (q/ (x)), while 
\P* is the pseudoscalar meson state with valence quarks of flavour / and /'. 
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that thanks to the choice r$ = —r^ in our MA setup the lattice axial current (747^75(73 is 
conserved (only broken by soft mass terms) and hence the matrix elements of the operator 

Qi in Dfl nglc 34 and ^Double 1234 vanish as ( M 34)2 _ ^ 3 + ^ in the Hmit ^ _^ Q 

This implies that no double pole occurs in Da. At non- vanishing lattice spacing there 
exists, however, an 0(a 2 ) single pole contribution in .D a mgle 12 because, owing to = 
r\, the lattice axial current i?27^75<7i is broken by discretization effects (see Ref. ifTTl and 
Appendix A of Ref. CD). 

For the case j 7^ 1, when no similar GB-pole simplifications can occur, double GB- 
pole terms strongly suppressed (like l/{p 2 ) 2 ) at large p 2 are to be expected in Dy(p). 
However, precisely owing to this strong suppression in practice, within our statistical errors 
and in the ranges of quark masses and p 2 we explore (see Table 2 of Ref. [47 ] and sect. IB.6I ). 
we hardly see in our lattice data any effects that can reliably be ascribed to double GB-pole 
contributions. On the contrary, we do find clear numerical evidence for single GB-pole 
contributions, which indeed at high p 2 are only suppressed as 1/p 2 . We thus decided to 
ignore double GB-pole terms in our valence mass chiral extrapolations. 

This choice is also justified a posteriori by the results of the valence chiral fits based on 
the ansatz (IB. 17b - A subset of these results is illustrated in Fig. [6] There we display typical 
examples of the effect of GB-pole subtractions in the matrix elements of the dynamics 
matrix at two values of /3. As can be seen, after the subtraction, a smooth dependence 
upon fj, va i (or equivalently on M 2 2 ) is observed. Combining the valence chiral limit lattice 
estimator of Dy and Z q , we are able to get reliable estimates of the intermediate quantities 
Zjf(p 2 ; a 2 p 2 ;0,a^ sea ). 

B.3 Sea chiral limit 

At fixed f3 and a 2 p 2 , we fit Zlj t (p 2 ;a 2 p 2 ;0,a/j lsea ) data to a first order polynomial in 
a 2 fi 2 ea . This choice is dictated by the expectation that effects of spontaneous chiral sym- 
metry, which may induce a dependence on |/i seo |, are strongly suppressed, and in practice 
immaterial within errors, in quantities like our RC-estimators that are evaluated at momen- 
tum scales p 2 S> Aq CD - In fact we find that the dependence on the sea quark mass is hardly 
visible within our statistical error bars, as shown in Fig. [7] Moreover, we have checked that 
repeating the whole analysis using a constant fit leads to similar RC results, though affected 
by smaller errors and often, but not always, yielding acceptable x 2 ' s - Hence, we conser- 
vatively decided to perform the sea chiral extrapolation using a linear fit in a 2 [i 2 ea . We 
construct in this way the RC estimators Zf^ (p 2 ; a 2 p 2 ; 0, 0). 
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Figure 6: GB-pole subtraction and valence chiral limit of D23, -D33, D44 and D55 plotted 
versus clival, for (3 = 3.9, a[x sea = 0.0040 and (ap) 2 m 1.565 (left column) and (3 = 
4.05, afi sea = 0.0030 and (ap) 2 rj 1.568 (right column). 
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Figure 7: The quantities Z^' (p 2 ; a 2 p 2 ;0, a^i sea ) and Z^' (p 2 ; a 2 p 2 ; 0, afi sea ), taken at the 
valence chiral limit, as functions of a 2 n 2 ea , for a typical lattice momentum choice (see 
inset) giving a 2 p 2 ~ 1.56, for the three (3 values considered in this paper. 



B.4 Removal of 0(a 2 g 2 ) cutoff effects 



We will obtain improved chiral limit RC estimators, Z^ 1 _impr (p 2 ; a 2 p 2 ), by removing from 
our Zjj 1 (p 2 ; a 2 p 2 ; 0, 0) lattice data perturbative discretization errors. This can be done up 
to 0(a 2 g 2 ) exploiting the one-loop perturbative results obtained ll53l l54l in the massless 
lattice theory for the quark propagator form factor Si, related to the quark-field RC by 
Z g (p) = Si(p), and the dynamics matrix elements, v.Lz. 
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-9 p2 



+ 0(a 4 g 2 ,9 4 ' 
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(B.19) 



The values of the coefficients Cq \ k = 1,2,3 can be found in Eq. (34) of Ref. B71 . while 



the values of the coefficients b\- and c)- are available in Refs. [53ll54l. 

In the numerical evaluation of the perturbative corrections, we take the coupling con- 
stant g 2 as the simple boosted coupling g 2 = g^/ (P) . For the average plaquette (P) 
we employ the non-perturbative values [0.5689,0.5825,0.6014] for /3 = [3.8,3.9,4.05], 
respectively. The important impact of the perturbative corrections in removing the un- 
wanted a 2 p 2 dependence is illustrated, for the case of /3 = 3.8, in Fig. [8j In this figure, 
the uncorrected values of Z^ 1 (//g = cl(/3)~ 2 ; a 2 p 2 ; 0, 0) are compared with the values of 
a(/3) -2 ; a 2 p 2 ) obtained setting either g 2 = g 2 , = 6//3 or (as we did in 



the end) g 2 = g 2 . 
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• Uncorrected 
■ 0(gg a 2 ) 

♦ 0(gV> 



1 

1 } 

5 5 



(a) 



Uncorrected 

Oltl a") 
0(gV) 



I 



(b) 



Figure 8: The effect of subtracting from Z^Vo = a(/3)" 2 ; a 2 p 2 ; 0, 0) and Z^'Oo = 
a((3)~ 2 ; a 2 p 2 ; 0, 0) at (3 = 3.8 (biue dots,) the 0(a 2 g 2 ) correction, setting either g 2 = 
(red squares) or g 2 = g 2 (green diamonds). 



B.5 Absence of wrong chirality mixings 

In Fig. |9j one can clearly see that for all the operators of interest the mixing coefficients 
Ay are very small (in fact vanishing within errors in the range of p 2 that we eventually 
use for extracting RCs). We also find that this is systematically more and more so as (3 
increases, well in line with our expectation that in our lattice setup wrong chirality mixing 
effects are reduced to mere 0(a 2 ) artifacts. For these reasons the effects of A have been 
neglected in our final RC analysis, where we have assumed a fully continuum-like relation 
between renormalized and bare operators. In addition we checked that repeating the whole 
analysis with the tiny effects of A on the relation (1B.1I) properly taken into account leads to 
no significant changes in the values of RCs. 

B.6 Final RC estimates from Ml and M2 method 

Having extrapolated the (improved) RCs estimators to the valence and sea chiral limit at 
each value of the momenta, we evolve Z^j 1 -impr (p 2 ; a 2 p 2 ) from the scale p 2 to a com- 
mon scale /jg by using the known matrix formula for the NLO running of the operators 
Qi ir5Tl[T5ll . obtaining Z? 1 _impr (/XQ; a 2 p 2 ). This step is necessary in order to disentangle 
the 0(a 2 p 2 ) cutoff effects from the genuine continuum p 2 dependence. Notice also that the 
actual value of fiQ has no impact on the RGI results of the RCs. As is customary, we take 
= a" 1 (/?) for each 0, with a" 1 (3.8, 3.9, 4.05) = [2.0, 2.3, 2.9] GeV. 
Of course we still allow for a residual dependence on a 2 p 2 . In order to deal with these 
cutoff effects, following Ref. R71l . we use two methods. Method Ml consists in fitting 
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Figure 9: The behaviour of the mixing coefficients Ajj, as a function of a?p 2 for (3 = 4.05. 



r^RI impr^2. a 2p2^ tQ jjjg \ mQW ansa tz 

z Ri'- impr(/i 2. a 2 f) = Z f-^^D + Xij (ap) 2 (B.20) 

in the large momentum region, 1.0 < a 2 p 2 < 2.2. As expected, the dependence of the 
slopes on f3 is very mild and not significant within our statistical errors. We thus treat 
Ajj as a /3-independent quantity. 

According to the ansatz (IB .201 ) a simultaneous linear extrapolation in a 2 p 2 = at each 
P was performed. The extrapolated values, zf"^ ~ impr (/iQ), are finally used to evaluate via 
the NLO running matrix formula of Ref. H3, the quantities Zff s (Ml) and Z^ GI (M1). 
Therefore, the MS scheme we use here is the one defined by Buras et al. in Ref. lTT5l . This 
definition of the MS scheme, which has become standard, differs from the one of Ref. II5TT1 
proposed by Ciuchini et al. in the treatment of the four-fermion evanescent operators ap- 
pearing in the calculation of the two-loop anomalous dimensions. 

In Fig. [10] the simultaneous best linear fits in a 2 p 2 at our three /3's of Zij are shown. 
We recall that both in the analysis and in the figures of this Appendix, only data points 
corresponding to the momenta p satisfying the constraint (IB.16I ) are used and shown. 

The idea of the M2 method is instead to separately average at each j3 the values of 
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Z^j 1 _impr (//o! a 2 p 2 ) over a narrow interval of momenta (ideally just one point), which 
has to be kept fixed in physical units for all /3's. We have chosen this interval to be p 2 € 
[8.0, 9.5 ] GeV 2 . In this way, at the price of giving up the reduction of cutoff effects implied 
by the Ml method, no assumptions are introduced in the RC analysis about the detailed 
form of lattice artifacts and/or the adequacy of NLO anomalous dimensions to describe the 
RC-evolution at scales below p 2 ~ 9 GeV 2 . 

The a 2 -scaling of renormalized quantities (in this work operator matrix elements) con- 
structed using RCs determined with the M2 method will of course be in general different 
from the one of their Ml method counterparts, but the continuum limit results for these 
quantities, if attainable from both methods with controlled errors, should be consistent with 
each other (see e.g. Appendix E). 
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Figure 10: Z™' ~ impT (fj% = a((3)~ 2 ; a 2 p 2 ) for {ij} = {22,23,32,33,44,45,54,55} as 
functions of a 2 p 2 for the three (3 values considered in our study. The straight lines represent 
the simultaneous linear fit to the lattice data in the interval 1.0 < a 2 p 2 < 2.2 at the three 
P's. 39 



C Renormalization Constant results 



In Tables [9] and \\0\ we collect values of Zp and Z$ calculated in the RI-MOM scheme. 
Results are obtained with methods Ml and M2 PTTl at each value of the gauge coupling in 
MS and RI-MOM at 2 GeV. We have used the three-loop conversion formula from RI-MOM 
to MS 1551. 



RC(M1) 


P = 3.80 


P = 3.90 


P = 4.05 


MS at 2 GeV 


Z P 


0.411(12) 


0.437(07) 


0.477(06) 


Z S 


0.734(09) 


0.713(05) 


0.699(06) 


RI-MOM at 2 GeV 


Z P 


0.338(10) 


0.359(06) 


0.391(05) 


Z s 


0.603(07) 


0.586(04) 


0.574(05) 



Table 9: Z P and Z s results, using the Ml method at P = 3.80, 3.90 and 4.05 in MS and 
RI-MOM at 2 GeV. 



RC(M2) 


P = 3.80 


P = 3.90 


P = 4.05 


MS at 2 GeV 


z P 


0.531(05) 


0.518(06) 


0.520(04) 


Z S 


0.813(07) 


0.776(06) 


0.735(04) 


RI-MOM at 2 GeV 


z P 


0.437(04) 


0.426(05) 


0.427(04) 


Z s 


0.668(06) 


0.637(05) 


0.603(04) 



Table 10: Z P and Z s results, using the M2 method, at P = 3.80, 3.90 and 4.05 in MS and 
RI-MOM at 2 GeV. 

The RC matrices of the four-fermion operators Zq, (c.f. Eq. (IA.10I )). are listed below. 
We present results obtained from both Ml and M2 methods, in MS and RI-MOM at 2 GeV. 
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(MS, 2 GeV): 



Z Q (/3 = 3.80;M1) 



Z Q (/3 = 3.90;M1) 



Zq(/3 = 4.05;M1) 



Z Q (/3 = 3.80;M2) 



Z Q (/3 = 3.90;M2) 



/ 0.421(13) 





/ 0.439(07) 




V 

/ 0.493(06) 




V o 

/ 0.440(09) 







/ 0.447(05) 







Z Q (/3 = 4.05;M2) 



/ 0.495(05) 









0.490(13 
0.021(02 







0.504(07 
0.023(01 







0.552(08 
0.025(01 







0.512(07 
0.043(02 







0.514(05 
0.038(01 







0.555(05 
0.033(01 









0.209(07) 
0.224(10) 







0.206(05) 
0.250(05) 







0.221(07) 
0.284(07) 







0.272(05) 
0.382(05) 







0.260(04) 
0.362(04) 







0.264(04) 
0.343(04) 











0.255(08) 
-0.245(09) 





0.280(05) 
-0.246(05) 





0.318(05) 
-0.261(08) 





0.367(05) 
-0.148(04) 





0.361(04) 
-0.169(03) 





0.360(03) 
-0.214(04) 







-0.006(02) 
0.600(15) 





-0.007(01) 
0.620(08) 





-0.004(01) 
0.690(11) 





-0.005(02) 
0.547(09) 





-0.005(01) 
0.576(05) 





-0.003(01) 
0.663(08) 
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(RI-MOM, 2 GeV): 



Z Q (/3 = 3.80;M1) 



Z Q (/3 = 3.90;M1) 



Zq(/3 = 4.05;M1) 



Z Q (/3 = 3.80;M2) 



Z Q (/3 = 3.90;M2) 



/ 0.415(12) 





/ 0.432(07) 




V 

/ 0.486(06) 




V o 

/ 0.433(08) 







/ 0.441(04) 







Z Q (/3 = 4.05;M2) 



/ 0.487(05) 









0.503(13 
0.016(01 







0.517(07 
0.018(01 







0.566(08 
0.019(01 







0.527(07 
0.034(01 







0.528(05 
0.031(01 







0.570(05 
0.026(01 









0.237(08) 
0.190(08) 







0.237(05) 
0.212(05) 







0.256(07) 
0.241(06) 







0.318(05) 
0.324(04) 







0.304(04) 
0.307(04) 







0.306(05) 
0.291(03) 











0.236(08) 
-0.239(08) 





0.259(05) 
-0.241(05) 





0.294(05) 
-0.256(07) 





0.338(04) 
-0.149(04) 





0.332(03) 
-0.169(03) 





0.331(03) 
-0.212(04) 







-0.013(02) 
0.572(14) 





-0.014(01) 
0.591(08) 





-0.012(01) 
0.659(10) 





-0.011(02) 
0.522(09) 





-0.012(01) 
0.550(05) 





-0.011(01) 
0.632(08) 



42 



D Lattice data on masses and matrix elements 



In the following Tables we gather our bare results at all values of /3 and combinations of 
quark masses for (i) pseudoscalar meson masses and pseudoscalar meson decay constants 
in lattice units (Tables [TT1 [T2l and [TBI already published in Ref. [8 ]; (ii) the ratio of the (bare) 
four-fermion operators = {K \Oi\K )/{K \Oi\K ) (i = 2, . . . , 5) (see Tables [11 [15] 
and Eg); (iii) the quantities £ (i = 2, . . . , 5) (see Tables El [l8]and[T9]|. 



(3 = 3.80 (24 3 x 48)a 4 



afl£ — a^i se a 


afi« s » 


aAf 34 


aM 12 


aF 34 


aF 12 




0.0165 


0.2558(08) 


0.3393(25) 


0.0894(04) 


0.0883(15) 


0.0080 


0.0200 


0.2731(07) 


0.3532(23) 


0.0913(04) 


0.0895(15) 




0.0250 


0.2961(07) 


0.3718(21) 


0.0936(04) 


0.0909(15) 




0.0165 


0.2712(04) 


0.3508(16) 


0.0924(03) 


0.0900(16) 


0.0110 


0.0200 


0.2877(04) 


0.3644(14) 


0.0942(03) 


0.0910(16) 




0.0250 


0.3098(04) 


0.3828(12) 


0.0966(03) 


0.0924(16) 



Table 11: Pseudoscalar masses and decay constants at f3 = 3.80. 



P = 4.05 (32 3 x 64)a 4 



am = afi sea 




aM 34 


aM 12 


aF 34 


aF 12 


0.0030 


0.0120 


0.1602(08) 


0.1931(18) 


0.0564(03) 


0.0558(07) 


0.0150 


0.1751(08) 


0.2053(17) 


0.0578(03) 


0.0566(07) 


0.0180 


0.1889(08) 


0.2169(16) 


0.0591(03) 


0.0573(07) 


0.0060 


0.0120 


0.1739(06) 


0.2034(11) 


0.0600(04) 


0.0585(08) 


0.0150 


0.1877(06) 


0.2153(11) 


0.0613(03) 


0.0596(08) 


0.0180 


0.2007(06) 


0.2266(11) 


0.0625(03) 


0.0605(08) 


0.0080 


0.0120 


0.1840(05) 


0.2127(09) 


0.0615(04) 


0.0604(12) 


0.0150 


0.1972(05) 


0.2242(09) 


0.0627(04) 


0.0616(12) 


0.0180 


0.2097(05) 


0.2351(09) 


0.0638(04) 


0.0626(12) 



Table 12: Pseudoscalar masses and decay constants at j3 = 4.05. 
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P = 3.90 (24 3 x 48)a 4 


afi£ = an sea 




aM 34 


aM 12 


aF 34 


Z7>12 

ar 


0.0040 


0.0150 


0.2060(05) 


0.2639(11) 


0.0724(03) 


0.0705(09) 


0.0220 


0.2401(05) 


0.2915(11) 


0.0757(03) 


0.0725(09) 


0.0270 


0.2619(05) 


0.3096(11) 


0.0777(03) 


0.0759(09) 


0.0064 


0.0150 


0.2179(08) 


0.2762(16) 


0.0755(05) 


0.0736(10) 


0.0220 


0.2506(07) 


0.3028(15) 


0.0785(05) 


0.0759(09) 


0.0270 


0.2717(07) 


0.3204(14) 


0.0805(04) 


0.0774(09) 


0.0085 


0.0150 


0.2283(07) 


0.2849(17) 


0.0773(03) 


0.0755(09) 


0.0220 


0.2598(07) 


0.3109(15) 


0.0804(03) 


0.0779(08) 


0.0270 


0.2803(07) 


0.3281(15) 


0.0823(03) 


0.0794(09) 


0.0100 


0.0150 


0.2351(07) 


0.2892(14) 


0.0784(04) 


0.0761(09) 


0.0220 


0.2659(07) 


0.3154(12) 


0.0815(04) 


0.0787(08) 


0.0270 


0.2860(06) 


0.3328(12) 


0.0834(04) 


0.0802(09) 


P = 3.90 (32 3 x 64)o 4 


0.0030 


0.0150 


0.1982(04) 


0.2558(13) 


0.0720(03) 


0.0701(07) 


0.0220 


0.2329(04) 


0.2838(11) 


0.0750(03) 


0.0720(08) 


0.0270 


0.2550(04) 


0.3021(11) 


0.0770(03) 


0.0745(08) 


0.0040 


0.0150 


0.2041(04) 


0.2644(15) 


0.0727(03) 


0.0702(11) 


0.0220 


0.2381(04) 


0.2917(17) 


0.0758(03) 


0.0722(10) 


0.0270 


0.2599(04) 


0.3096(15) 


0.0777(03) 


0.0753(09) 



Table 13: Pseudoscalar masses and decay constants at P = 3.90. 
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P = 3.80 (24 3 x 48)a 4 



CL/J,£ = a^ sea 






4" 


„(*>) 

it 4 


*?> 




0.0165 


-13.14(07) 


3.19(02) 


24.29(12) 


8.18(04) 


0.0080 


0.0200 


-11.92(06) 


2.89(01) 


21.93(10) 


7.43(03) 




0.0250 


-10.58(04) 


2.56(01) 


19.34(08) 


6.61(03) 




0.0165 


-12.16(06) 


2.95(01) 


21.89(10) 


7.41(03) 


0.0110 


0.0200 


-11.15(05) 


2.70(01) 


19.99(08) 


6.81(03) 




0.0250 


-9.99(04) 


2.41(01) 


17.84(07) 


6.13(02) 



Table 14: R- for i = 2, . . . , 5, as obtained from Eq. (14.171 ). at each combination of the 
quark mass pair (a/i^, a/i« s >') an d at /3 = 3.80. 



(3 = 4.05 (32 3 x 64)a 4 



afi£ — a^ sea 














0.0120 


-20.47(12) 


5.07(03) 


40.49(22) 


13.32(07) 


0.0030 


0.0150 


-17.35(09) 


4.29(02) 


34.29(17) 


11.36(06) 




0.0180 


-15.10(07) 


3.72(02) 


29.78(13) 


9.93(04) 




0.0120 


-16.61(09) 


4.11(02) 


32.49(15) 


10.79(05) 


0.0060 


0.0150 


-14.55(07) 


3.59(02) 


28.35(11) 


9.48(04) 




0.0180 


-12.97(05) 


3.20(02) 


25.19(09) 


8.48(03) 




0.0120 


-15.13(05) 


3.74(01) 


29.36(09) 


9.79(03) 


0.0080 


0.0150 


-13.42(04) 


3.31(01) 


25.99(07) 


8.72(02) 




0.0180 


-12.09(03) 


2.97(01) 


23.34(06) 


7.88(02) 



Table 15: R- for i = 2, . . . , 5, as obtained from Eq. (14. 17b . at each combination of the 
quark mass pair (afi£, afi« s » ) and at (3 = 4.05. 
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P = 3.90 (24 3 x 48)a 4 


afl£ = a^ sea 




K 2 




K 4 




0.0040 


0.0150 


-16.42(06) 


4.15(01) 


32.51(12) 


10.79(04) 


0.0220 


-13.11(04) 


3.20(01) 


24.99(07) 


8.41(02) 


0.0270 


-11.38(03) 


2.77(01) 


21.57(06) 


7.32(02) 


0.0064 


0.0150 


-15.37(08) 


3.77(02) 


29.02(14) 


9.68(04) 


0.0220 


-12.21(05) 


2.98(01) 


22.89(09) 


7.75(03) 


0.0270 


-10.70(04) 


2.60(01) 


19.97(07) 


6.82(02) 


0.0085 


0.0150 


-14.10(05) 


3.44(01) 


26.56(09) 


8.92(03) 


0.0220 


-11.41(03) 


2.77(01) 


21.32(06) 


7.25(02) 


0.0270 


-10.09(03) 


2.45(01) 


18.75(05) 


6.44(01) 


0.0100 


0.0150 


-13.52(05) 


3.31(01) 


25.30(09) 


8.53(03) 


0.0220 


-11.07(04) 


2.70(01) 


20.55(06) 


7.02(02) 


0.0270 


-9.85(03) 


2.39(01) 


18.28(05) 


6.26(02) 


P = 3.90 (32 3 x 64)a 4 


0.0030 


0.0150 


-17.18(08) 


4.22(02) 


33.45(14) 


11.09(05) 


0.0220 


-13.27(05) 


3.24(01) 


25.66(12) 


8.63(03) 


0.0270 


-11.47(05) 


2.80(01) 


22.09(09) 


7.50(03) 


0.0040 


0.0150 


-16.32(09) 


4.10(02) 


32.22(18) 


10.68(06) 


0.0220 


-13.01(06) 


3.17(01) 


24.88(12) 


8.34(04) 


0.0270 


-11.31(05) 


2.74(01) 


21.42(09) 


7.27(03) 



Table 16: K- for i = 2, . . . , 5, as obtained from Eq. (14.171 ). at each combination of the 
quark mass pair (afie, a/i« s ") and at /3 = 3.90. 
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j3 = 3.80 (24 3 x 48)a 4 



afi£ — afi sea 










tA 0) 




0.0165 


1.015(12) 


0.247(03) 


1.877(21) 


0.632(07) 


0.0080 


0.0200 


1.029(11) 


0.249(02) 


1.892(21) 


0.641(07) 




0.0250 


1.046(11) 


0.253(02) 


1.912(20) 


0.654(07) 




0.0165 


1.038(07) 


0.252(02) 


1.868(12) 


0.632(04) 


0.0110 


0.0200 


1.050(07) 


0.254(02) 


1.883(11) 


0.641(04) 




0.0250 


1.065(06) 


0.257(02) 


1.902(11) 


0.654(04) 



Table 17: Bare ' for i = 2, . . . , 5, as obtained from Eq. (I4.15I ). at each combination of 
the quark mass pair (a/j,g, a[i" s ») and at /3 = 3.80. 



P = 4.05 (32 3 x 64)a 4 



afj,£ = a/i sea 














0.0120 


0.915(10) 


0.227(02) 


1.810(17) 


0.596(06) 


0.0030 


0.0150 


0.929(09) 


0.230(02) 


1.837(17) 


0.609(06) 




0.0180 


0.943(09) 


0.232(02) 


1.860(16) 


0.620(05) 




0.0120 


0.931(09) 


0.230(02) 


1.820(18) 


0.605(06) 


0.0060 


0.0150 


0.946(09) 


0.234(02) 


1.843(17) 


0.617(06) 




0.0180 


0.960(09) 


0.237(02) 


1.864(17) 


0.628(06) 




0.0120 


0.948(09) 


0.234(02) 


1.842(17) 


0.614(06) 


0.0080 


0.0150 


0.962(08) 


0.237(02) 


1.861(16) 


0.625(05) 




0.0180 


0.974(08) 


0.240(02) 


1.8883(16) 


0.636(05) 



Table 18: Bare £j_B- for i = 2, . . . , 5, as obtained from Eq. (14- 15b . at each combination of 
the quark mass pair (a/j,g, a[i« s ") and at /3 = 4.05. 
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P = 3.90 (24 3 x 48)a 4 


a\n = ap, sea 












0.0040 


0.0150 


0.961(07) 


0.236(02) 


1.848(14) 


0.613(05) 


0.0220 


0.991(07) 


0.242(02) 


1.888(13) 


0.635(04) 


0.0270 


1.009(07) 


0.245(02) 


1.911(13) 


0.649(04) 


0.0064 


0.0150 


0.979(09) 


0.240(02) 


1.848(20) 


0.617(06) 


0.0220 


1.006(09) 


0.245(02) 


1.887(16) 


0.639(05) 


0.0270 


1.023(08) 


0.249(02) 


1.909(16) 


0.652(05) 


0.0085 


0.0150 


0.987(08) 


0.241(02) 


1.860(16) 


0.624(05) 


0.0220 


1.014(08) 


0.246(02) 


1.894(15) 


0.644(05) 


0.0270 


1.030(08) 


0.250(02) 


1.914(15) 


0.657(05) 


0.0100 


0.0150 


0.993(12) 


0.243(03) 


1.861(24) 


0.627(08) 


0.0220 


1.019(12) 


0.248(03) 


1.894(23) 


0.647(08) 


0.0270 


1.036(12) 


0.252(03) 


1.915(23) 


0.660(08) 


P = 3.90 (32 3 x 64)a 4 


0.0030 


0.0150 


0.953(05) 


0.234(01) 


1.855(10) 


0.615(03) 


0.0220 


0.982(05) 


0.240(01) 


1.898(09) 


0.638(03) 


0.0270 


1.000(05) 


0.244(01) 


1.924(09) 


0.653(03) 


0.0040 


0.0150 


0.962(06) 


0.236(02) 


1.831(12) 


0.608(04) 


0.0220 


0.992(06) 


0.242(02) 


1.875(12) 


0.632(04) 


0.0270 


1.011(06) 


0.246(02) 


1.901(11) 


0.646(04) 



Table 19: Bare ^iB\ ' for i = 2, . . . , 5, as obtained from Eq. (14.151 ). at each combination of 
the quark mass pair (a/ie, a/i" s ") and at f3 = 3.90. 
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E Results for R4 and Bi 



In this appendix we present in detail our results in the MS scheme of Ref. [15] at 2 GeV 
for the quantities Ri and Bi (c.f. Eqs. (14.191 14.201 ) and Eq. (14.181 ) respectively). We also 
give the Ri results computed in the indirect way of Eq. (I5.6I ). In Table 1201 we gather results 
obtained employing Ml-type RCs and using ChPT (NLO) fit formula, polynomial and linear 
fit functions (see Eqs. (I5.2H5.3I ). and n = 2 and n = 1 of Eq. (15.11 ) respectively). In Table [2T1 
we show the respective results when using M2-type RCs. Instead of using the definition of 
Eq. ( 14.191 ), we have employed a slightly different but equivalent one which reads 

5/ _ 5 [ G k G k]\mi m . . 

l ~ 1 [G%G%]\ M2 {hA) 



where indices Ml and M2 refer to the use of the respective type of renormalisation constants 
and we define G K ' \(mi,M2) = (A/i, A/2)- We find that the quantity de- 

fined in Eq. (IE. II ) has smaller 0(a 2 ) effects. 

In Figs. [TT] and [12] we show the combined fit for the ratios, R\ and bag parameters Bi 
(i = 2, . . . , 5) against the light quark mass when M2-type RCs are used. 

We remark that a good agreement between the continuum limit results for the bag 
parameters Bi and the matrix elements ratios Ri obtained using M2-type RCs and their 
counterparts based on Ml-type RCs, as we observe for i = 1,2,3,4 (for i = 1 see also 
Ref. [8]]), provides a valuable check of the smallness of residual systematic errors in the 
evaluation of RI-MOM RCs with the Ml -method. In particular possible systematic errors 
stemming in the Ml -method from the inadequacy at non-high momenta (p 2 ) of the perturba- 
tive operator anomalous dimensions used in the analysis or from the removal of the leading 
cutoff effects via a linear fit in p 2 are strongly reduced or absent when using the M2-method 
for RCs. This is so because in this latter approach (see Ref. [47] and Appendix [B]) the RCs 
are extracted from Landau gauge correlators at a rather high p 2 -value (fixed to ~ 9 GeV 2 
for all /3's) but comes at the price of generically larger lattice artifacts on the RCs, which 
we partly suppress by removing the perturbatively known 0(a 2 g 2 ) contributions. For the 
case of Oi with i = 1, 2, 3, 4 the resulting cutoff effects on Ri and Bi (see Figs ITTI and [P21 
appear to be under control and the continuum extrapolation is reliable. On the contrary, the 
case of B§ and R^ (see panel d) of the figures above) is a typical one where too large cutoff 
effects affecting the M2-type RCs make unreliable the results appearing in the i = 5-lines 
of Table Ell 
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Fit 


i 


Ri 


i?i(via Eq. lETS])) 


Bi 


ChPT 


1 


1 


1 


0.51(02) 


2 


-14.3(04) 


-16.2(24) 


0.53(02) 


3 


5.1(03) 


5.7(09) 


0.93(07) 


4 


26.3(09) 


29.8(45) 


0.83(04) 


5 


7.0(04) 


7.7(12) 


0.64(06) 




Pol 


1 


1 


1 


0.52(02) 


2 


-15.0(07) 


-16.5(26) 


0.54(04) 


3 


5.3(04) 


5.8(11) 


0.95(12) 


4 


26.0(14) 


30.0(48) 


0.82(06) 


5 


6.9(07) 


7.8(16) 


0.63(09) 




L 


1 


1 


1 


0.53(01) 


2 


-14.3(04) 


-16.5(24) 


0.54(02) 


3 


5.1(03) 


5.8(09) 


0.95(07) 


4 


24.9(08) 


29.1(44) 


0.81(04) 


5 


6.6(04) 


7.5(12) 


0.62(06) 



Table 20: Ri (direct computation through Eq. ( 14.191 ) and indirect computation through 
Eq. (15.61 )) and B t results using Ml-type RCs for three kinds of fit function, namely a ChPT 
(NLO) fit, a polynomial and a linear fit with respect to the light quark mass. For % = 2, 3 
the ChPT (NLO) fit formula for Ri coincides with the linear one (we refer to results of the 
3rd column). 
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Fit 


i 


Ri 


Ri(via Eq. (J5HD) 


Bi 


ChPT 


1 


1 


1 


0.52(02) 


2 


-14.5(04) 


-15.7(22) 


0.51(01) 


3 


4.9(02) 


5.4(08) 


0.88(03) 


4 


27.1(08) 


28.4(41) 


0.79(02) 


5 


5.0(02) 


5.9(09) 


0.49(03) 




Pol 


1 


1 


1 


0.52(02) 


2 


-15.4(07) 


-16.1(24) 


0.53(02) 


3 


5.1(03) 


5.5(08) 


0.90(05) 


4 


27.4(14) 


28.6(43) 


0.78(04) 


5 


4.8(05) 


5.9(10) 


0.48(05) 




L 


1 


1 


1 


0.53(01) 


2 


-14.5(04) 


-16.0(23) 


0.52(01) 


3 


4.9(02) 


5.5(08) 


0.90(03) 


4 


25.6(07) 


28.0(40) 


0.78(02) 


5 


4.7(03) 


5.8(09) 


0.48(03) 



Table 21: Ri (direct computation through Eq. (IE. II ) and indirect computation through 
Eq. (15.61 )) and B t results using M2-type RCs for three kinds of fit function, namely a ChPT 
(NLO) fit, a polynomial and a linear fit with respect to the light quark mass. For i = 2, 3 
the ChPT (NLO) fit formula for Ri coincides with the linear one (we refer to results of the 
3rd column). The results in the lines corresponding to i = 5 here are not reliable, due to 
very large cutoff effects resulting in this case from the use of M2-type RCs (see text). 
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Finally we give our continuum results for Bi and Ri in the MS scheme of Buras et al., 
defined in Ref. 03, and the RI-MOM scheme at 3 GeV, see Tables l22l and Irrespectively. 



MS (3 GeV) 



Bi 


B 2 


B 3 




B 5 


0.51(02) 


0.51(02) 


0.85(07) 


0.82(04) 


0.66(07) 




Ri 


R2 


R 3 




R 5 


1 


-16.3(06) 


5.5(04) 


30.6(13) 


8.2(05) 



Table 22: Continuum limit results for Bi and Ri, renormalized in the MS scheme of 
Ref. HI at 3 GeV. 



RI-MOM (3 GeV) 



Bi 


B 2 


Bi 


B A 


B 5 


0.50(02) 


0.63(03) 


1.07(09) 


0.95(06) 


0.75(09) 




Ri 


R2 


Ri 




Rs 


1 


-15.4(06) 


5.3(03) 


26.9(12) 


7.1(05) 



Table 23: Continuum limit results for Bi and Ri, renormalized in the RI-MOM scheme at 
3 GeV. 
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0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
(a) 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
(c) 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 

rofii 

(b) 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 
(d) 



Figure 11: Solid lines in panels (a) and (b) show the behaviour vs. the renormalized 
light quark mass of the combined chiral and continuum fits (according to the polynomial 
formula if5.il) with n = 2) of the i?- (see Eq. (E.l\) ). with i = 2 and i = 3 respectively, 
renormalized in the MS scheme of Ref. f73l7 at 2 GeV with the M2-type RCs. In panels (c) 
and (d), solid lines, instead, show the combined chiral and continuum described by NLO- 
ChPT, Eq. &5.3\) for i = 4 and i = 5, respectively. The full black line is the continuum 
limit curve. The dashed black line represents the continuum limit curve in the case of the 
linear fit ansatz. Black open circles and triangles stand for the results at the physical point 
corresponding to the polynomial (panels (a) and (b)) and ChPT fit (panels (c) and (d)), and 
linear fit ansatz, respectively. 
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1.20 - 
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1 1.00 - 
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0.70 
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P = 3.90 
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ChPT Fit 
Lin. Fit 
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0. 
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° 
S 
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0. 
0. 
0. 
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40 
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20 
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P = 3.90 
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-5- 



-l- 



-i — t- 



-f — -r 

_L i- 



0.02 0.04 0.06 0.08 0.10 0.12 0.14 

rah 

(d) 



Figure 12: Solid lines in panels (a) to (d) show the behaviour vs. the renormalized light 
quark mass of the combined chiral and continuum fits (according to the ChPT fit formula 
(5.2\) ) for the Bi parameters with i = 2, . . . , 5 respectively, renormalized in the MS scheme 
of Ref. ifT5l? at 2 GeV with the M2-type RCs. The full black line is the continuum limit 
curve (f5. JD . The dashed black line represents the continuum limit curve in the case of the 
linear fit ansatz. Black open circles and triangles stand for the results at the physical point 
corresponding to the ChPT fit and linear fit ansatz, respectively. 
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